Elastic energy loss with respect to the reaction plane in a Monte-Carlo model 
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We present a computation of tt° nuclear modification factor with respect to the reaction plane in 
Au+Au collisions at ^/sjvjv = 200 GeV, based on a Monte-Carlo model of elastic energy loss of hard 
partons traversing the bulk hydrodynamical medium created in ultrarelativistic heavy-ion collisions. 
We find the incoherent nature of elastic energy loss incompatible with the measured data. 
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I. INTRODUCTION 

The interaction of hard high transverse momentum (pt) 
partons with the bulk medium created in heavy-ion colli- 
sions, leading to partonic energy loss and jet quenching, 
is regarded as one of the most important tools to study 
the properties of the medium. Initially, medium-induced 
gluon radiation [U-IH has been regarded as the most im- 
portant energy loss mechanism. However, it was pointed 
out that also elastic energy loss (where the energy from 
a leading parton is carried away in the recoil of a scat- 
tering partner) is a sizeable contribution in perturbativc 
Quantum Chromodynamics (pQCD) prl— Till] - In particu- 
lar, elastic energy loss seems to help explaining the data 
on high Pt single electron suppression, widely believed to 
be driven by heavy-quark energy loss [l2l . [l3j , for which 
the radiative energy loss is small due to the dead cone 
effect [lH ■ This raises the question of the importance of 
elastic processes to light-quark energy loss. 
Consistently, in microscopic calculations, a large elastic 
contribution to the total energy loss was found 
under the assumption that the degrees of freedom (DOF) 
in the medium are perturbatively interacting partons 
with a thermal spectrum or (almost) free thermal quasi- 
particles. In contrast, in a phenomenological top-down 
analysis of the pathlength dependence of the energy loss 
seen in the measured data, based on the fact that in 
a constant medium elastic energy loss increases linearly 
with pathlength whereas radiative energy loss increases 
quadratically, a contribution of more than 10% elastic 
energy loss to the total was ruled out [18J. These re- 
sults, however, are not in manifest contradiction, as the 
assumption of quasi-free scattering partners being avail- 
able in the medium may not be realized in nature, and 
as until recently [l6| no comparison with strongly path- 
length dependent observables have been performed in mi- 
croscopical models of elastic energy loss. 
The aim of this paper is to resolve this discrepancy by ap- 
plying a detailed microscopical Monte Carlo (MC) model 



of elastic energy loss, presented in Ref. [17|, to the nu- 
clear suppression factor as a function of the angle of 
outgoing hadron relative to the reaction plane, -R^A (</>), 
for different collision centralities. Our strategy is as fol- 
lows: We adjust the parameters of our MC energy loss 
model such that for a given hydrodynamical description, 
Raa{4>) data is reproduced for 0-10% central collisions. 
Without adjusting further parameters, we then apply the 
model to non-central collisions. Going from central to pe- 
ripheral collisions, the mean density of the system drops 
as well as the mean pathlength a parton travels in the 
medium. Both effects decrease the amount of suppres- 
sion (or increase the suppression factor Raa)> but a lin- 
ear weighting of the pathlength should lead to less effect 
than a quadratic weighting. Thus, in an elastic energy 
loss model, or any incoherent mechanism, the mean Raa 
increases more slowly from central to peripheral colli- 
sions than in a coherent radiative energy loss model. At 
the same time, for a given centrality the split between 
emission in the reaction plane and out of the reaction 
plane is a function of weighted pathlength, so we expect 
a larger spread for radiative energy loss than for elas- 
tic energy loss. Thus, by comparison with the PHENIX 
data 19, 20], we expect to be able to probe whether the 
characteristic scaling of elastic energy loss is realized. 
The manuscript is organized as follows: First, we describe 
the hydrodynamical model used for the evolution of the 
bulk matter. We then give the essential outline of our MC 
model for elastic parton-medium interactions, followed by 
a presentation of the results and a discussion. 



II. HYDRODYNAMICS 

We describe the bulk matter as a locally thermalized 
fluid. For this purpose, we solve the ideal fluid hydro- 
dynamic equations 



dpT» v = 0, 



(1) 
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where T^ v = (e + P)u^u v — g^ v P is the stress-energy 
tensor and e is energy density, P is pressure and is 
the fluid flow four-velocity. We also need an Equation 
of State (EoS) which relates pressure with the energy 
density and net-baryon number density, P = P(e, tlb)- 
Our choice here is the EoS from Laine and Schroder 12111 . 
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Since we consider here particle production only at mid- 
rapidity, we assume the net-baryon density to be negligi- 
ble. For the same reason, we assume longitudinal boost- 
invariance. We solve the (2+l)-dimensional numerical 
problem using the SHASTA algorithm [H|23j]. 
As an initial state we use the sWN profile from Ref . [3] ■ 
This profile reproduces the centrality dependence of 
dN/dr] at mid-rapidity 25] 1 . Overall normalization is 
fixed so that we get the same initial entropy in central 
collisions as in Ref. E3l- Motivated by the EKRT model 

0.17 fm. 



[28j and Ref. [27J , we set the initial time to To 
Centrality classesjire defined using the optical Glauber 
model as in Ref. 



27] 



Thermal spectra for hadrons are calculated using the con- 
ventional Cooper- Frye method [29j , where particle emis- 
sion is calculated from a constant-temperature surface. 
The freeze-out temperature Td ec = 160 MeV is fixed so 
that we reproduce the measured px spectrum of pions 
30] up to p T ss 2 GeV at various centralities [25| . Strong 
and electromagnetic two- and three-particle decays are 
taken into account before comparing with experimental 
data. 

Elliptic flow, which measures the momentum anisotropy 
of final state particles, is produced during hydrodynami- 
cal evolution since the initial state has spatial anisotropy. 
Because of this anisotropy, also the average pathlength of 
hard partons is angle dependent. So, to make reasonable 
energy loss calculations, it is important to constrain the 
initial eccentricity from the bulk matter elliptic flow. In 
Fig. [T] we show the integrated elliptic flow. We do not get 
enough elliptic flow in central collisions, but, as explained 
in [32j , this is due to the lack of initial state density fluc- 
tuations and a different reference plane definition than 
in the experimental analysis. These results show that we 
can reproduce the bulk observables reasonably well, so 
that we can do meaningful energy loss calculations using 
this hydrodynamical model as background. 
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FIG. 1: (Color online) Integrated elliptic flow for Au+Au 
collisions at -y/sjvjv = 200 GeV. Data are from the PHOBOS 
collaboration [3l| . The data are shown with statistical and 
systematic errors added in quadrature. 



u, d, s, u, d, s, g in the initial state, and over all possible 
parton type pairs (kl) in the final state. In general, the 
scattering rate depends on the frame, and in particular 
on the high-energy parton's 4-momentum pi, on the flow 
4-velocity u(x) and on the temperature T(x) of the fluid 
at each space-time location x. 

In the local rest-frame of the fluid, we can express the 
scattering rate as follows [TtJ : 



III. THE MONTE CARLO SIMULATION 

We model the energy loss of a hard parton by incoherent 
partonic 2—^2 processes in pQCD, with scattering part- 
ners sampled from the medium. Our simulation of en- 
ergy losses of high-energy partons in the produced QCD 
matter is based on the scattering rate for a high-energy 
parton of a type i, 

r i (p 1 ,u(x),T(x)) = ^2 Fij^ki(pi,u(x),T(x)), (2) 

j(kl) 

where we account for all possible partonic processes ij — ¥ 
kl by summing over all types of collision partners, j = 



When comparing the computed -»= — ' 7o with the data, the same 
Glauber model must be applied, see Ref. |2(J . 



Fij-ykl = ^„ 2 g2 / 2 dEfafj {E2, T)f2jj_>fc2 {E\ , E2, TU 2 ), 

(3) 

where 

fl l:i ^ki(Ei,E2,m 2 ) = dslsvij^kiis)}- (4) 

J2m? 

Here E\ is the energy of the high-energy parton i in 
this frame and E 2 is the energy of the thermal parti- 
cle j with a distribution function fj [E2 , T) , which is 
the Bose-Einstein distribution f g = fB(E 2 ,T) for glu- 
ons and the Fermi-Dirac distribution f q = fjj(E2,T) for 
quarks. The scattering cross section aij^ki (s) appearing 
in (j4j depends on the standard Mandelstam variable s. 
A thermal-mass-like overall cut-off scale m = s m g s T is 
introduced in order to regularize the singularities appear- 
ing in the cross section when the momentum exchange 
between partons approaches zero. Here g s is the strong 
coupling constant and s m is a parameter of the order of 
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one. Since we cannot meaningfully control the running 

2 

of a s = we keep the strong coupling constant fixed 
with momentum scale. 

To initiate the hard massless parton of a type i in each 
event, we sample the LO pQCD single-jet production 
spectrum (for more details, see (33[) at yi = and 
a t PTmin Pt < \fs/2. For the parton distribution 
functions (PDFs), we use the CTEQ6L1 set [11. For 
now, we neglect the nuclear effects to the PDFs |35l - l37| . 
since these are small in comparison with the ones aris- 
ing from the final state interactions with the medium. 
The initial rapidity yi is randomly generated in the range 
[ymin, ymax] from a flat distribution. This fixes the hard- 
parton energy E and polar angle 9 of its momentum vec- 
tor p = (p x ,py,p z ). The azimuth angle 4>, defined with 
the reaction plane, is evenly distributed between [0,27r]. 
The hard parton is assumed to start interacting with the 
medium at the initial longitudinal proper time t of our 
hydrodynamical model. Since in the cm. frame all hard 
partons are produced in the Lorentz-contracted overlap 
region at z 0, the longitudinal position at later times 
(before the first collision at r > To) is assumed to be de- 
termined by the longitudinal momentum only. The initial 
time and longitudinal coordinates for the hard parton are 
thus to = To cosh yi and zq = To sinh yi . The coordinates 
on the transverse plane in the beginning of the simulation 
are then xq = Xi + Pg-to and yo = yi + j^t , where the 
parton position in the transverse plane at t = 0, (xi,yi), 
is sampled from the nuclear overlap function 



T AA (h) = d 2 sT A (s + h/2)T A ( S -b/2), 



(5) 



where b is the impact parameter. The nuclear thickness 
function T A (s) is defined as usual, 



Ta (s) = J dzp A (s,z) 



(0) 



where for the nuclear density p A {v) we use the standard 

Woods-Saxon distribution for the gold nucleus. 

The hard parton propagates through the plasma in small 

time steps At, during which we propagate the parton 

in position space. The probability for not colliding in 

this time interval is assumed to be given by the Poisson 

distribution 



the partial scattering rates (|3]) to the total scattering rate 
@. After scattering, the final state parton with highest 
energy is chosen as the new hard parton to be propagated 
further, for which we repeat the procedure outlined above 
with the next timestep. 

We take into account the system's slow transformation 
from quark-gluon plasma to hadron gas by using an ef- 
fective temperature 



rp ( 30 

T e ff = 2 € 



1/4 



(8) 



where gQ ~ g g + ^2Nfg q = 4p is the quark-gluon plasma 
degrees of freedom with gluon and quark DOF being g g = 
16 and g q = 6, respectively, and number of quark flavors 
Nf = 3. Our simulation ends at T — 165 GeV, when the 
system can be considered to have completely transformed 
into hadron gas. We assume no significant interaction 
between the high-energy parton and the fully hadronic 
medium. 

The outcome of the procedure described above is a 
medium-modified distribution of high-energy partons, 
' x In order to calculate the nuclear modifica- 



dp T dy 

tion factor for neutral pions, 



R AA (P T ,y,b) 



dNl A /dP T dy 



T AA (b)daPP/dP T dy' 



(9) 



we hadronize this distribution by convoluting the ob- 
tained partonic distribution with the fragmentation func- 
tion Df^. n (z,/J,p): 



dN 



AA ^TT°+X 



dP T dy 



E 



dp T dy 



dN AA ^f+ x 
dp T dy 



I dzDf^. n (z, n 2 F )6(P T - zp T ), 
Jo 



(10) 



where z — — is the fraction of the final parton mo- 
mentum px available to the hadron with momentum Pt , 
and fj,F ~ Pt is the fragmentation scale. The Kniehl- 
Kramer-Potter (KKP) fragmentation functions [38[ have 
been used for this work. 



P(No collisions in At) = e 



-TiAt 



(7) 



where Ti is the total scattering rate ((2]) for the hard 
parton of the type i. Hence the probability to col- 
lide at least once during the time At is 1 — e~ FiAt « 
TiAt + 0((TiAt) 2 ). For small enough At we can as- 
sume that there will be at most one collision. Note that 
since we calculate the scattering rates (HJ) in the local rest 
frame of the quark-gluon plasma fluid element, we must 
boost also the time step At to the same frame. Should 
a scattering happen, the probability Pij^ki for a given 
type of scattering process is determined by the ratios of 



IV. RESULTS 

In the following, our interest is in the neutral pions with 
high transverse momentum PT min > 5.0 GeV, to account 
for the observation that at low Pt the main hadron pro- 
duction mechanism is not independent fragmentation. 
We average all observables across the rapidity window 
[—0.35,0.35] which corresponds to the PHENIX accep- 
tance. Due to the non-eikonal propagation of the hard 
partons (see discussion in Ref. [17]) in our simulation, 
the initial rapidity window is y m in = —ymax = — 1 in 
order to account for all the possible partons falling into 
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FIG. 2: (Color online) Left panel: The tv° nuclear modification factor, averaged over the reaction plane angle, for 0-10% 
centrality. Middle and right panel: The it nuclear modification factor dependence on the reaction plane angle Acj> for 40-50% 
(middle panel) and 50-60% centrality (right panel). Only the smallest- and the largest-angle bins are shown. The strong 
coupling constant is a s = 0.5 and mass parameter s m = 1.0. The simulation points (solid squares and triangles) are connected 
with lines to guide the eye. PHENIX data is from [20J (0-10% centrality, open squares) and [19( (40-50% and 50-60% centrality, 
open triangles). Colored bars with small cap represent statistical errors; black bars with wide cap are systematic errors. 



the final rapidity window. As an attempt to emulate the 
reaction plane measurements by PHENIX [19], the high- 
momentum partons within the final rapidity acceptance 
have been divided in six bins of angle Acf> with respect to 
the reaction plane, from A(j> = - 15° to A0 = 75 - 90°. 



For this work, the value of strong coupling constant is 
chosen to be a s = 0.5, which produces roughly the right 
amount of suppression to match the measured nuclear 
modification in the 0-10% centrality bin. In effect, such 
tuning emulates also the incoherent higher-order pro- 
cesses, whose contribution can be significant. The idea 
is then to examine how well our tuned MC model agrees 
with the measurement in other centralities. As our inter- 
est here is particularly in the dependence on the reaction 
plane angle, we concentrate on the more peripheral cases 
where the effect should be visible best. The simulation 
results for 0-10%, 40-50% and 50-60% centrality bins, 
compared with the measurements, are shown in Fig. [5J 



It is clear from the figure that our model cannot repro- 
duce the reaction plane angle dependence seen in the 
PHENIX experiment. Also the inclusive, angle-averaged 
nuclear modification factor fails to match with the exper- 
imental data: The computed suppression decreases too 
slowly as one advances to the more peripheral collisions. 



V. DISCUSSION 

Our result, obtained for Raa{<P) hi a microscopical 
well-controlled MC model with a realistic hydro back- 
ground, is an important one since it demonstrates that a 
purely incoherent energy-loss framework seems to be in 
a striking contradiction with the present RHIC data. In 
essence, we confirm the findings of the phenomenological 
analysis [HI that the pathlength dependence of elastic 
energy loss is incompatible with the data. Taken to- 
gether with evidence that even the quadratic pathlength 
dependence of coherent radiative energy loss in pQCD 
may face some challenge in describing the data [39|], this 
implies that any contribution of elastic energy loss must 
be rather small (we refrain from trying a combined fit of 
coherent radiative and incoherent elastic energy loss to 
the data at this point, but 10%, as obtained in Ref. [l8| . 
seems certainly reasonable). 

The implication is that the main assumption made in 
pQCD calculations of elastic energy loss, i.e. that the 
medium DOF are almost free (quasi)-particles which can 
take a sizeable amount of recoil energy away from a lead- 
ing parton may not be true in nature. The assumption 
of static color dipoles (leading to purely radiative energy 
loss), difficult as it may be to justify, seems to work much 
better in explaining the data. Given that also a cubic 
pathlength dependence is consistent with the data [ioT |. 
it is unfortunately too early to draw a firm conclusion 
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with regard to the DOF and the energy loss mechanism 
realized in the medium; more systematic studies need to 
be done. However, given that an ideal fluid description 
of the bulk medium, which has a zero mean free path 
assumption implicitly made, works well to describe the 
bulk, it should perhaps not come as a surprise that an 
infinite mean free path assumption (the ideal parton gas 
approximation for fj in Eq. does not work well to 
describe the same medium as seen by a hard parton. 
At this point, it is also not clear what the implications of 
our findings for the question of heavy-quark energy loss 
are. The reason is that heavy quarks are rare probes, i.e. 
they are not thermally excited in the medium. While 
a reaction like qg — > qg in which the incoming quark is 
hard, but the outgoing gluon carries 90% of its energy 
must be counted as 10% energy loss from the leading 
parton if the incoming quark is light, a c quark is always 
uniquely identified by its flavour and the same reaction 
must be interpreted as 90% energy loss from the tagged 
heavy quark. Thus, elastic energy loss is generically more 
efficient for heavy quarks than for light quarks, and we 
will investigate the question of heavy quark elastic energy 
loss in a future work. 



In this work the initial state density fluctuations were 
neglected. Because of this, all the results were deter- 
mined using the reaction plane, which is defined by the 
impact parameter. However, the experimental results for 
Raa^) are calculated with respect to event plane, which 
may differ from the reaction plane. At least at low pt 
this difference is important when calculating vi [32| . To 
study whether the initial state density fluctuations and a 
reference plane definition consistent with the experiment 
are as important in Raa(<P) is l e ft f° r future work. 
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